On the possible sources of gravitational wave bursts detectable today 
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We discuss the possibility that galactic gravitational wave sources might give burst signals at 
a rate of several events per year, detectable by state-of-the-art detectors. We are stimulated by 
the results of the data collected by the EXPLORER and NAUTILUS bar detectors in the 2001 
run, which suggest an excess of coincidences between the two detectors, when the resonant bars 
are orthogonal to the galactic plane. Signals due to the coalescence of galactic compact binaries 
fulfill the energy requirements but are problematic for lack of known candidates with the necessary 
merging rate. We examine the limits imposed by galactic dynamics on the mass loss of the Galaxy 
due to GW emission, and we use them to put constraints also on the GW radiation from exotic 
objects, like binaries made of primordial black holes. We discuss the possibility that the events 
are due to GW bursts coming repeatedly from a single or a few compact sources. We examine 
different possible realizations of this idea, such as accreting neutron stars, strange quark stars, and 
the highly magnetized neutron stars ( "magnetars") introduced to explain Soft Gamma Repeaters. 
Various possibilities are excluded or appear very unlikely, while others at present cannot be excluded. 
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I. INTRODUCTION 

In ref. [1] the ROG collaboration has presented the 
analysis of the data collected by the EXPLORER and 
NAUTILUS resonant bars during nine months in the year 
2001. When the number of coincidences between the two 
resonant bars is plotted against sidereal time, one finds 
an excess of events with respect to the expected back- 
ground, concentrated around sidereal hour four. At this 
sidereal hour the two bars are oriented perpendicularly 
to the galactic plane, and therefore their sensitivity for 
galactic sources of gravitational waves (GWs) is maxi- 
mal. Furthermore, for these events the energy deposed in 
one bar is well correlated with the energy deposed in the 
other. The significance of this observation has been de- 
bated [2-4] , and certainly further experimental work will 
be necessary to put the indications on a firmer ground. 
New data from bars and interferometers, from long data 
taking runs, should soon clarify the situation. 

As we will discuss below, from a theoretical point 
of view the existence of GW bursts with the in- 
tensity and the rate necessary to explain the EX- 
PLORER/NAUTILUS results, if they will be confirmed 
by further data, would certainly be a great surprise. It 
would then be necessary to reconsider with an open mind 
a number of unexpected possibilities. The purpose of this 
paper is to provide a general framework for such a study, 
indicating what are the difficulties, which directions of 
investigation appear to be more promising and which are 
ruled out. In particular, we hope that our work will help 
to set the stage for the analysis of future data. 

The structure of the paper is as follows. In sect. II we 
present some of the basic aspects to be explained, that is: 
(1) the energy which must be released in GWs in order to 
explain the events observed, and (2) the rate which is in- 
ferred from the EXPLORER/NAUTILUS observations. 
The other crucial experimental information is the distri- 
bution of the events as a function of sidereal time. This is 



presented in sect. Ill, together with a detailed discussion 
of what can be learned from sidereal time analysis. We 
think that this section has also a more general method- 
ological interest. In sect. IV we examine one of the most 
obvious candidates, the coalescence of a compact binary 
system. We find that compact objects of solar masses, 
at typical galactic distances, would account for the en- 
ergy release. However, the rate expected from binaries 
made of neutron stars and/or black holes is many order 
of magnitude smaller than the observed rate. We then 
turn our attention to more exotic binary systems, such 
as those involving primordial black holes. The theoreti- 
cal uncertainties on exotic events make more difficult to 
put direct limits on their rate. However, the coalescence 
of binary systems should be a phenomenon which takes 
place at a more or less steady rate for a time comparable 
to the age of the Galaxy. We show in sect. V that the 
corresponding loss of mass of the Galaxy into GWs is 
constrained by galactic dynamics. Again, we think that 
this section is of more general interest, independently of 
the application to the EXPLORER/NAUTILUS data. 

These considerations lead us to suggest, in sect. VI, 
that the all events might be generated by a single (or 
just a few) "GW burster", i.e., by an object which emits 
repeatedly bursts of GWs, and we examine different pos- 
sible implementations of this idea. 

In sect. VII we show that the signal cannot be ac- 
counted for by non-gravitational phenomena like the de- 
position of recoil energy in the bars due to the pas- 
sage of massive particles. We present our conclusions 
in sect. VIII. 



II. ENERGY REQUIREMENTS 

Given a GW with Fourier spectra h+ lX (f), the gravi- 
tational energy radiated per unit area and unit frequency 
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is given by [5] [88] 
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For a wave coming from an arbitrary direction and with 
arbitrary polarization, a detector does not measure di- 
rectly ft+(/) and h x (f) but rather the combination 
h(f) = F+h+(f) + F x h x (f), where F +;X arc the de- 
tector pattern function. For a bar, F + = sin 2 9 cos 2ip, 
F x = sin 2 6 sin 2tp, where ip describes the polarizations 
and 6 is the angle of arrival, measured from the bar axis. 
The events we consider here are recorded when the bar 
is orthogonal to the galactic plane, and therefore when 
sin 2 6 ~ 1. Assuming that the source is emitting ran- 
domly polarized GWs, we average over ip. Then 
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where (. . .) denotes the average over the polariza- 
tions and we used (sin 2 2-0) = (cos 2 2-0) = 1/2 and 
(sin 2ip cos 2ip) = 0. Therefore eq. (1) becomes 
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For a burst, we assume that (ft 2 (/)) is equal to a con- 
stant ft, 2 between a frequency f m - m and a frequency / max . 
Denoting by r the distance to the source and assum- 
ing isotropic emission, the total radiated energy A£ ra( j 
is therefore given by 

47r 2 r 2 c 3^2 j 
AE rad = - (/max - /min) • ( 4 ) 

Typically, the term / min is negligible in comparison with 
/max) an d therefore we write simply 
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For EXPLORER and NAUTILUS the value of h c is re- 
lated to the energy E s deposed in the bar by 



h r = 2.5 x Kr 21 Hz -1 
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Therefore, under the hypothesis that the events recorded 
indeed correspond to GWs, each burst originates from a 
process that liberated in GWs the energy 
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A typical value of E s for the considered events is E s <~ 
100 mK for the 2001 data, see ref. [1]. 

In terms of the GW amplitude ft, again assuming a flat 
Fourier spectrum up to a frequency / max , for a burst of 



duration At ~ l// max we have ft ~ h c f maK , so eq. (6) 
gives 



ft ~ 2.5 x 10" 
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Assuming that a couple of events are due to backgrounds, 
as is expected in a two hours period [1] , we estimate that 
the 8 events recorded correspond to a signal rate of order 
200 events/yr. The challenge is therefore to explain what 
kind of source could give such a strong GW emission, 
at such a high rate. We also remark that there is no 
observed neutrino counterpart of these events [6] . 



III. SIDEREAL TIME ANALYSIS 

The other crucial piece of information is the distribu- 
tion of the number of coincidences when plotted against 
sidereal time, as stressed in particular in refs. [7-9] . The 
purpose of this section is to discuss in detail what we can 
learn about the location (and possibly the polarization) 
of the sources from the sidereal time analysis. 

The experimental data for the 2001 run are shown in 
fig. 1, adapted from ref. [1]. The coincidences are binned 
in one hour bins corresponding to the sidereal hour of 
arrival. The number of accidental coincidences (dotted 
line) is estimated shifting the data stream of one detector 
with respect to the data stream of the other by a step 
St = 2 s and measuring the number of coincidences, which 
now are all accidentals. The analysis is repeated with a 
step 2St, then with 3St, etc., up to 100<5t. If n(j) is the 
number of coincidences found with a shift jSt, then the 
average number of accidental coincidences is taken to be 
n = (1/100) n(j). We see from the figure that there 

is an excess of coincidences with respect to the expected 
number of accidentals, at sidereal hours 3 and 4. The 
data in fig. 1 refers to events collected only in periods 
with more than 12 hr of continuous operation. If we 
consider all periods with more than 1 hr of continuous 
operation the events at the peak becomes 4 + 4 rather 
that 4 + 3 as in fig. 1, and the background also rises 
slightly, see fig. 7 of ref. [1]. 

While the statistical significance of the peak is not 
large, two important facts make this result more intrigu- 
ing: (1) sidereal hour 4 is a very special moment. In fact, 
a peak at this value of sidereal time is predicted to appear 
for sources in the galactic plane. This is a very general 
prediction, independent on the precise location and na- 
ture of the sources. As we will discuss in detail below, 
depending on the distribution of sources in the galactic 
plane (e.g., a uniform distribution in the disk, sources 
concentrated in the galactic center, etc.) a second peak, 
possibly smaller, can appear at another value of sidereal 
time, but the value of t of this second peak depends in- 
stead strongly on the specific position of the sources. (2) 
For the 8 events at the peak, the energy deposed in one 
bar is very well correlated with the energy deposed in the 



3 



5_ 
4_ 

3_ 
2_. 



1 



FIG. 1: The number of coincidences found in the 2001 
run (solid line) and the estimated background (dashed line), 
against sidereal time i = (£ex P i + iNau)/2 (adapted from fig. 5 
of ref. [1]). 

other, while for the events at other sidereal hours the en- 
ergy deposed in one bar and in the other are completely 
uncorrelated [1]. 

To extract physical informations from the plot against 
sidereal hour we proceed as follows. Let 6 denote the an- 
gle between the direction of a source and the longitudinal 
axis of the bar. The response of a single bar to GWs from 
this source is 

h(f) = sin 2 8 (h+ cos 2ip + h x sin 2ipj . (9) 

Therefore the response in amplitude is proportional to 
sin 2 6 and, since the energy is quadratic in the amplitude, 
the response of the detector in energy is proportional to 
sin 4 8. Furthermore, there can be an effect due to the 
polarization of the waves, i.e., to the angle ip. We will 
consider first the case of unpolarized GWs and then we 
will discuss the possible modification due to the depen- 
dence on ip. 

A. Randomly polarized GWs 

If the GWs come from an ensemble of sources or from a 
single source which emits randomly polarized waves, the 
effect of ip will be averaged out and the energy deposed 
in a detector will just be proportional to sin 4 8. Of course 
6 = 6(t) because of the rotation of the Earth. Using a bit 
of geometry (whose details are left to appendix A) one 
can see that 

cos 9{t) — \n x cos(a — fit) + n y sin(a — Qt)] cos 5 + 

+n z sin<5, (10) 

where a, S are the right ascension and declination of the 
source, a is the azimuth of the bar, I its latitude, n x = 
— cos a sin I, n y = sin a, n z — cos a cos/, ft = 2-7r/24 is 
the rotation frequency of the Earth in units 1 / (sidereal 
hours), t is the local time measured in sidereal hours, and 
all angles, including a, S, are expressed in radians. 

The axis of the two detectors are aligned to within 
a few degrees of one another, so that the response of 
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FIG. 2: The location of the galactic plane in the equatorial co- 
ordinate system. The position of the galactic center is marked 
by a circle. 



the two bars to any source can be considered identical. 
The latitudes and azimuths arc: for EXPLORER, I = 
46.457V and a = 39° E while for NAUTILUS I = 41.827V, 
a = M°E. The longitudes are instead 6.20E and 12.67E, 
respectively. Given their difference in longitude, the local 
sidereal time at the EXPLORER and at the NAUTILUS 
locations differ by At ~ 25.88 min. We will denote by 
t the average between the two local sidereal times, t = 
(£ex P i + ^Nau)/2 — £ex P i + 0.216 hr. (The same variable 
was used in ref. [1]. To compare with other experiments 
it might be more convenient to use Greenwich sidereal 
time ^Green, related to i by t ~ icroon + 0.629 hr). 

In the equatorial coordinate system (a, S) the galactic 
plane is represented by the curve in fig. 2. For sources 
in the galactic plane a plot of sin 4 9 against sidereal time 
t, with < t < 24, always shows two maxima: one peak 
is close to sidereal hour 4, because at that moment the 
bar turns out to be almost perpendicular to the galac- 
tic plane. The precise position of the source within the 
galactic plane results only in a minor variation of the pre- 
cise value of t ~ 4 where sin 2 6 = 1. Instead, the sidereal 
time at which there is the second peak depends strongly 
on the location of the source. 

Fig. 3 shows sin 4 6 as a function of t for a source 
located in the galactic center, i.e., a = 266.405° and 
5 = —28.936°. In fig. 4 we show sin 4 6 for another source 
in the galactic plane, taking as an example a source 
with equatorial coordinates (a = 135°, 5 = —40°). Both 
curves have a maximum in the bin 4 < t < 5 (more pre- 
cisely, at i ~ 4.55 for the galactic center and t ~ 4.27 
for the source of fig. 4) and they have a second peak 
whose position depends on the source location. Fur- 
thermore, they have a minimum where sin 4 6 becomes 
zero. Moving away from the galactic plane, e.g. increas- 
ing the declination S at fixed a (see fig. 2) at first the 
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FIG. 3: sin 4 9 against sidereal time t (in hr) for the galactic 
center. 




FIG. 4: sin 4 6 against sidereal time t (in hr) for a source with 
coordinates a = 135° and 8 — —40°. 



curve still has two maxima, with positions that depend 
on the values of (a, 6), but then one arrives at a criti- 
cal value of 5 where the two maxima coalesce. Above 
this critical value there is just a single maximum, and 
as we increase further S the value at the peak becomes 
smaller than one. At the same time the value at the min- 
imum moves away from zero, and therefore the response 
curve becomes much flatter. Fig. 5 shows the curve for a 
source with a = 266.405°, S = 80°. In the limit 6 = 90° 
we see from eq. (10) that the curve becomes fiat, at a 
value which depends on a, I, and which in our case is 
sin 4 ~ 0.504. A plot of sin 2 as a function of the lo- 
cal sidereal time t and of the declination S is shown in 
fig. 6. The effect of changing a on this figure is simply 
to produce a shift in the origin of sidereal time, as we see 
from the fact that in eq. (10) a and t appear only in the 
combination a — fit. 

We now discuss how these response curves are reflected 
in a plot of the number of observed coincidences, N c , 
against sidereal time. Consider for definiteness the sit- 
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FIG. 5: sin 4 6 against sidereal time t (in hr) for a source with 
a = 266.405°, <5 = 80°. 



FIG. 6: sin 2 6 as a function of sidereal time t and of the 
declination 5, for a = 135°. 



uation in which the sources are in the direction of the 
galactic center. In fig. 3 we have shown sin 4 as a func- 
tion of i for these sources, and it is intuitive that the 
maxima of N c , as a function of i, coincide with the two 
maxima of sin 4 0, since at that values of sidereal time 
the bars have their best orientation with respect to these 
sources. However, it is important to understand that the 
functional form of N c , and in particular the width of the 
maxima, is not the same as that of the function sin 4 0. 
Rather, it depends crucially on the ratio between the 
energy of the events and the energy threshold of the de- 
tectors [4, 7-9]. To understand this point, suppose that 
in a given observation time arrive n events in each one- 
hour bin and restrict for the moment to the simplified 
situation in which a signal arriving at a generic sidereal 
time t deposes in the bar an energy E sin 4 0(t), with E 
the same for all signals. Consider first the limiting case 
in which the threshold of the detector is very low with 
respect to Eq, e.g., is represented by the dotted line A in 
fig. 7. In this case all the events are above threshold and 
therefore are recorded, except those falling exactly in the 
blind direction of the detectors, which in this example 
corresponds to the bins between sidereal hours 19 and 
23. Therefore in the bins corresponding to sidereal hours 
19 to 23 we would record zero events, while in all others 
we would observe n events per bin. In this limiting case, 
therefore, N c (t) shows no peak at all; rather, it is a flat 
function, N c (t) — n, except for a dip in correspondence 
to the blind direction. 

The opposite limiting case is represented by an energy 
threshold very close to E (the dotted line B in fig. 7) . In 
this case, only when we are in the two bins where sin 4 
becomes equal to one we can detect the signals, while 
as soon as we move to the next bin we fall below the 
threshold and no signal is detected. In this case N c = n 
in the bins 4 < t < 5 and 13 < t < 14, while N c — 
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FIG. 7: The energy of the events in the various sidereal time 
bins, compared to two possible thresholds of the detectors 
(the dotted lines marked as A and B). 



FIG. 8: A simulation of the number of events per bin, for 
a distribution of sources in the galactic plane (chosen to be 
uniform in galactic longitude) when the detector threshold is 
such that E/E = 0.9. 



otherwise. A plot of N c (t) will show two very narrow 
peaks, each one concentrated just in a single sidereal hour 
bin. 

In a more realistic situation, the above picture will be 
complicated by the fact that the signals will arrive with a 
distribution in energy, rather than being monoenergetic. 
Moreover, the energy of the events is in any case "dis- 
persed" by the fact that the events observed are a com- 
bination of the GW signal and of the noise. The above 
example, however, suffices to illustrate that the width 
of the peak is crucially affected by the energy distribu- 
tion of the events and by the detector thresholds, and 
can be very flat (for great values of the typical signal-to- 
noise ratio) to very narrow (for small SNR). One should 
also observe that, when the average SNR of the events 
is small, we can see sin 4 9 as a measure of the probabil- 
ity that, in a single detector, a signal is not lost below 
the energy threshold. In a two-detector correlation the 
probability that both detectors see the signal is therefore 
rather measured by sin 8 9, which therefore is a more ap- 
propriate envelope for the graph in fig. 7, and this results 
in an even narrower peak. 

It is also interesting to note that, in the case of large 
SNR, informations on the location of the sources cannot 
be extracted from N c (t), which is now basically flat, but 
it could be obtained plotting against sidereal time the av- 
erage energy of the coincident events, although probably 
a large sample of data would be needed, to compensate 
for the spread in the intrinsic energies of the signals. 

If we assume a given distribution of sources and a given 
typical value of the detector threshold with respect to the 
signal, we can simulate the distribution of the number of 
events as a function of sidereal time. In fig. 8 we show an 
example where, for illustration, we considered a popula- 
tion of sources all in the galactic plane, and distributed 
uniformly in galactic longitude (this is natural for sources 
seen from the Earth, if we can detect only sources which 
are not too far, say r < 0.5 kpc. If one would be sen- 



sible also to very far sources, then clearly toward the 
galactic center there would be many more sources than 
in the opposite direction). We assume that each source 
deposit in the bar an energy E$ sin 4 9(t), with Eq the 
same for all bursts and for all sources. We fix a detector 
threshold (in fig. 8 we have chosen a rather high thresh- 
old, E/Eq — 0.9). For each source we produce a plot as 
fig. 7, and we see which bins are above threshold. Then 
as the contribution of this source we take N c (t) = 1 for 
these bins and N c (t) — otherwise. Summing the contri- 
butions of all sources and normalizing the distribution so 
that the value at the peak is one, we get fig. 8. Of course 
this is a rather simplified model, but it illustrates that a 
distribution of galactic sources may produce a single peak 
close to sidereal hour four (whose width is controlled by 
the threshold that we set). 

Comparing with the experimental data in fig. I we 
see that the preliminary indication, within the statis- 
tics available, is that there could indeed be a peak near 
sidereal hour 4. No second peak of some statistical sig- 
nificance can be seen. From the above discussion, this 
points toward the existence of several sources distributed 
across the galactic disk. Independently of the details of 
the models considered in fig. 8, the result follows more 
generally from the fact that for a distribution of sources 
in the galactic plane each single source would contribute 
either to the bins around sidereal hour 4 common to all 
sources in the galactic plane, or to another bin which de- 
pends on the specific source location. Only in the former 
case the contributions from different sources add up and, 
within the present statistics, give rise to something which 
emerges over the background. 

However, here we should not rush toward conclusions. 
Even assuming that the events really correspond to GWs, 
it is clear that the full form of the curve can be under- 
stood only when sufficient statistics will be available. We 
will therefore keep open, for the moment, the possibility 
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that further structures in the sidereal time plot might 
appear with further data, and we will examine different 
possibilities. 

We also mention that the sidereal time analysis can be 
very useful even in the study of extragalactic sources [7, 
8]. However, eq. (7) clearly excludes an interpretation 
in terms of GWs of extragalactic origin for the EX- 
PLORER/NAUTILUS data, and therefore we will focus 
on galactic sources. 



B. Polarized GWs 

In sect. VI we will examine the possibility that all 
events come from just one source which emits repeat- 
edly GW bursts. One can expect from such a source a 
coherent motion of matter with some given quadrupolar 
pattern related to the source geometry and to its rotation 
axes, and consequently the emission of polarized GWs. 
As we will see in this section, the polarization can affect 
the sidereal time analysis in an interesting way. 

The response of the detector (in amplitude) is given 
in eq. (9), and the polarization angle tp depends on side- 
real time, just as 9. The calculation of the dependence of 
this response function on local sidereal time t, for a given 
location of the source and a given choice of the polariza- 
tion axes, is in principle straightforward but somewhat 
lengthy, and we give the details in Appendix A. 

In fig. 9 we show the result for a source located in 
the direction of the galactic center (a similar result is 
reported in ref. [10]). In this figure the envelope is sin 4 9, 
i.e. the response function (in energy) for an unpolarized 
source, and is therefore the same quantity already shown 
in fig. 3. The thick line is the function sin 4 9 cos 2 2ip, 
i.e. the response function (in energy) for a source which 
has only the + polarization, while the dotted curve is 
sin 4 #sin 2 2ip, i.e. the response function in energy for the 
x polarization. (The angle ip is measured with respect 
to an axis orthogonal to the propagation direction and 
lying along the galactic plane, see Appendix A). 

From this figure we see a number of interesting effects. 
First of all, the maxima are shifted with respect to the 
unpolarized case. This is due to the fact that the maxima 
of sin 2 2ij) and of cos 2 2ip do not coincide in general with 
the maxima of sin 4 9, and therefore the position of the 
peaks in the response function is determined by a combi- 
nation of the two effects. For the same reason, the values 
of the response function at the peaks will be smaller than 
one, and need not be equal among the two peaks. We also 
see that the peak close to sidereal hour 4 in the + po- 
larization is narrower compared to the unpolarized case, 
because for the + polarization the maximum of cos 2 20 
is close to the maximum of sin 4 9 (as we see from the 
fact that the position of the peak close to sidereal hour 
4 is only slightly shifted) and therefore the additional 
factor cos 2 2\p in the energy response gives a further sup- 
pression outside the peak. New blind directions appear, 
when cos 2 20 = or sin 2 2-0 = 0, which contribute to 




FIG. 9: The response function in energy for unpolarized GWs 
(the envelope), for the + polarization (thick line), and for 
the x polarization (dotted line), for a source at the Galactic 
center. 




FIG. 10: The square of the quantities shown in fig. 9 (which 
is the appropriate response function for a two-detector corre- 
lation) . 



make the peaks narrower. 

Recalling that the response function of a two-detector 
correlation, for low typical SNR, is the square of the 
single-detector response function, we see that finally we 
can obtain very narrow peaks, as shown in fig. 10. This 
means that in the case of a source emitting waves with 
a high degree of polarization along the galactic plane it 
is not necessary to have a very low value of the typical 
SNR in order to find a narrow peak in the number of ob- 
served coincidences N c (t). A threshold of order 0.5 in the 
response function of fig. 10 suffices to have N c (t) concen- 



7 




FIG. 11: The same as fig. 9, for a source with coordinates 
a = 135°, 5 = -40°, and with a shift if) -> V> - tpo with 
V>o = 0.25 rad. 



trated in just two hour bins around i = 4, for a source 
radiating the + polarization. Therefore the number of 
observed events is not necessarily a large underestimate 
of the total number of events, as it would instead be the 
case if most of the signals were lost below the threshold. 

In this example we have considered the (rather special) 
case in which the + polarization is perpendicular to the 
galactic plane, i.e., the two axes used to define h + and 
h x are one parallel and the other perpendicular to the 
galactic plane (see Appendix A). More in general, we can 
rotate these axes in the transverse plane by an angle tp . 
A source which, with respect to these new axes, emits 
purely the h+ (or purely the h x ) polarization will have 
a response function obtained shifting ip — + ip — tpo, and 
therefore the position of the peaks depends on tpo . As an 
example, we show in fig. 11 the response to the + and x 
polarization, for a source located at a = 135°, S = —40° 
(i.e., the same source as in fig. 4) and with ipQ = 0.25 
rad. 



IV. COALESCENCE OF COMPACT BINARIES 

As we will recall below, for binaries with neutron 
stars (NS) and/or black holes (BH) the expected galactic 
merging rates are far too small, compared to the value 
0(200) events per year discussed in sect. II. However, 
stimulated by the experimental data, one might wish 
to consider more unusual possibilities, like BHs of pri- 
mordial origin or other exotic objects. For this rea- 
son, and also in order to understand better the diffi- 
culties of finding a good candidate source for the EX- 
PLORER/NAUTILUS data, we start our discussion with 
the coalescence of compact binaries. 

The merging of NS-NS, NS-BH or BH-BH binaries is a 
process that generates GWs at the kHz, where the bars 
are sensitive. Furthermore, they typically liberate an en- 



ergy of order 10~ 2 M Q c 2 in GWs. This can be estimated 
first of all analytically, using the wave form of the inspiral 
phase, 

Mt) = -v (GMc)5/3 ^ (r) ( i± f^) C0S[$(T)] ' 

2 4/3 

M*) = -(GM c ) 5/3 ^ 3 (r) costsin[$(r)], (11) 

where M c is the chirp mass, M c = (M 1 M 2 ) 3/5 /(M 1 + 
M 2 ) 1/5 , M 1 ,M 2 are the masses of the two objects, i is 
the inclination of the orbit with respect to the line of 
sight, r is the time to coalescence, $ is the accumulated 
phase and ui gvJ (T) is the chirping frequency of the GW. 
Explicitly, 

^(t) = | (^pj (GM c )- 5 /V/ 8 + $ , (12) 



(r) 



( (GM c )- 5 / 8 r- 3 / 8 . 



v 4 



(13) 



Taking the Fourier transform in the saddle point approx- 
imation (see e.g. ref. [11]) one gets, apart from irrelevant 
phases, 



M/) 



1 / 5 \ 1/2 1 (GM C ) 5 / 6 fl + cos 2 t 



7T 2 / 3 V 24 



/ 7/e 



1 / 5 \ 1/2 1 (GM C ) 5 / 6 

Using eq. (1) and performing the angular integration one 
finds the energy spectrum, 



(15) 



Integrating up to a maximum frequency / max for which 
we are still in the inspiral phase we can estimate the total 
energy radiated during the inspiral phase, 

A^ad^^GM^/^ (16) 



or, inserting the numerical values, 
A£ rad ~ 4.2 x 10" 2 M Q c 2 



M \ 5/3 / f \ 2/3 

lv± r \ I J ma 



1.22M, 



1kHz 



(17) 

where the normalization of the chirp mass corresponds 
to M x = M 2 = 1.4M . 

This simple analytic treatment gives a value in very 
good agreement with that found numerically, from so- 
phisticated hydrodynamical simulations of NS-NS coales- 
cence including both the inspiral and plunge phases [12], 



A£ rad ~ (1-3) x 10~ 2 Mqc 2 



(18) 
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FIG. 12: The distribution of distances r against mass M for 
the eight 2001 events at the peak at sidereal hours 3 and 4. 
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FIG. 13: The distribution of distances r against cost for the 
eight 2001 events at the peak at sidereal hours 3 and 4, setting 
M = 1.35M©. 

Comparison with eq. (7) indicates that, from the en- 
ergetic point of view, the coalescence of compact bina- 
ries is an interesting candidate. Actually, eq. (7) has 
been obtained assuming, for lack of better informations, 
a square waveform in frequency space. Since in the case 
of an inspiral we have the exact waveform, we can be 
more precise. Averaging over the angle ip we have in fact 

^ = lsin 4 0(|M/)| 2 + IM/)| 2 ) , (19) 
while from eq. (14) we have 

ife4/)i 2 +ifex(/)i 2 = 24 j /3r2 (G ^r /3 ^), (20) 

where 

f 1 + cos 2 l\ 2 9 . . 

g{t)= ^ J +cos 2 t . (21) 

We can then substitute this expression for h c into eq. (6). 
We set / ~ 920 Hz (corresponding to the resonant fre- 
quency of the bars) . For events coming from the galactic 
plane, detected around sidereal hour four, we can also set 



sin 9 ~ 1 . Then we find 

E s ( M c \ 5/3 /8kpc\ 2 (gm 

lOOmK ' \1.22MqJ \ r ) V - 8 / 

(22) 

We have chosen as a reference value for g(t) its average 
value over the solid angle. Each event, i.e. each value of 
E s , determines a curve in the plane (M c ,r), where 

£lc = Mc (as) ^ (23) 

and this allows us to check whether each event can be 
interpreted with reasonable values for M c and r. The 
function g{i) varies between between 1/4 for cos i = 0, i.e. 
when the observer is in the plane of the orbit, and g{t) = 
2 for cos/, = 1, i.e. when the line of sight of the observer 
is perpendicular to plane of the orbit. Correspondingly, 
M c ranges between 0.50M C and 1.73M C . 

The results are shown in fig. 12, for the eight events at 
sidereal hours 3 and 4 in the 2001 data. We have written 
M c = M/2 1 / 5 , which is the chirp mass of a system with 
two equal masses Mi — M 2 = M, and we have fixed 
g{b) = 0.8. In fig. 13 we fix instead M = 1.35M Q , a 
typical value for a NS, and we plot r against cos i for the 
same eight events. 

We see that very natural values of M and r are ob- 
tained. Binary systems composed of compact solar-mass 
objects, at typical galactic distances, would be compati- 
ble with the energetic observed in the events. 

The problem with these sources, however, is the rate. 
Before the recent discovery of a new NS-NS binary, the 
rate of NS-NS coalescences in the Galaxy was estimated 
to be in the range 10~ 6 to 5x 10~ 4 mergings per year [13]. 
This estimate depends strongly on the shortest-lived sys- 
tems known, and the recent discovery of a new NS-NS 
binary, with the shortest known merging time (85 Myr), 
brings this estimate up by one order of magnitude, pos- 
sibly up to a factor 30 [14]; still, we are very far from 
the O(200) events per year needed to explain the NAU- 
TILUS/EXPLORER observations. The estimates of the 
coalescence rate of binary systems involving black holes 
(BH-BH and BH-NS) can be based only on stellar evolu- 
tion models, rather than on observations, since no such 
systems have yet been observed. The theoretical calcu- 
lations suffers from large uncertainties, but the rate for 
BH-BH and BH-NS coalescence is estimated to be of the 
same order of magnitude (or smaller) than the NS-NS 
coalescence rate, see e.g. ref. [15] and references therein. 

Under the stimulus of the NAUTILUS/EXPLORER 
data, one might consider more exotic possibilities, as for 
example black hole of primordial rather than astrophys- 
ical origin. In this case the estimates based on stellar 
evolution do not apply, and it is interesting to observe 
that primordial BHs produced in the early universe at 
the QCD phase transition would have today a mass of 
approximately 0.5M Q [16]. The density of MACHO's, 
and therefore also of primordial BHs of mass <~ 0.5M Q , 
can be constrained from the microlensing of stars from 
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the Large Magellanic Cloud [17, 18]. These data suggest 
that the fraction of our Galaxy's stellar mass which is in 
the form of black holes can be significantly larger than 
1%. With some assumptions, the fraction of primordial 
black holes in binary systems, and their distributions in 
eccentricity and semi-major axis has been estimated in 
[19, 20], and it has been found that the coalescence rate of 
primordial black hole binaries could be of order 5 x 1CP 2 
events per year per galaxy. This is still far too small 
to explain the data. Furthermore, a stellar population 
of primordial origin is expected to reside in the galac- 
tic halo rather than in the galactic disk, so there should 
be no reason to get an enhanced signal when the bars 
are oriented favorably with respect to the galactic plane. 
However, when one enters into the realm of exotic pos- 
sibilities there are large uncertainties in the theoretical 
estimates. For instance, for primordial black holes, there 
are large uncertainties on the detailed formation mech- 
anism, and therefore one should keep an open mind on 
this option. 

Other interesting informations might come observing 
that the coalescence of any population of compact bina- 
ries should be a phenomenon which takes place at a more 
or less steady rate for a timescale comparable to the age 
of the Galaxy. It is therefore useful to examine whether, 
independently of the nature of the sources, the corre- 
sponding mass loss of the Galaxy to GWs is compatible 
with known facts of galactic dynamics. We discuss this 
issue in the next section. 



V. UPPER LIMITS ON THE GW EMISSION 
FROM THE GALAXY 

From the energy radiated in each burst, eq. (7), and a 
rate 0(200) events/yr, we see that the mass lost by the 
Galaxy into GWs is of order 

-M~2^° f-^-V feV, (24) 
yr V8kpcy V lkHz / 

where now r is an average distance of these sources. If the 
events are indeed due to the merging of compact objects 
we know the energy which has been released by each 
burst, independently of the values of r and of / max , see 
cq. (18), and therefore we can write 

-M~(2-6)^. (25) 
yr 

This is a very large mass loss. For comparison, the 
mass loss of the Galaxy due to electromagnetic radia- 
tion is 9 x 10~ 3 M Q /yr. In this section we will investigate 
whether such a high value is compatible with known facts 
of galactic dynamics. 

If the mass loss is really due to GWs from a number 
of different sources distributed across the galactic disk, it 
is difficult to imagine that its rate changed dramatically 
with time over the age of the Galaxy. Assuming a steady 



mass loss at the level of eq. (25) for a time comparable 
to the age of the oldest clusters in the Galaxy, T <~ 1.2 x 
10 10 yr, gives a total mass loss over the history of the 
Galaxy of at least AM ~2x 10 10 M Q . For comparison, 
the total mass of the disk is estimated to be 6 x 10 10 Af Q 
and the total mass of the disk plus bulge and spheroid is 
9 x 10 10 M Q . The total halo mass is M ha i Q - 2 x 10 12 M Q 
[21], but for us this is a less relevant reference value, since 
this is the mass in a spherical halo with radius of order 
170 kpc, and therefore events coming from the halo would 
not show any special correlation with the galactic plane. 

Of course, the present value of the mass of the Galaxy 
gives only a first scale for comparison, since on the one 
hand we can imagine that the mass of the Galaxy could 
have been larger at the formation and, on the other hand, 
tighter constraints on the possible mass loss may come 
from galactic dynamics. Some of these issues have been 
considered many years ago [22-24], and it is interest- 
ing to go back to these considerations using the present 
knowledge of galactic dynamics, and to compare with the 
energy loss given in eq. (25). The most stringent limits 
are discussed in the following subsections. 

A. Effect of the mass loss on the radial velocity of 

stars 

Let us first recall the basic fact of galactic dynamics 
that, because of the differential rotation of the Galaxy, 
the radial velocities v r of stars as seen from the sun, at 
first order in the distance R between the sun and the star 
(after correcting for the solar motion and averaging over 
the peculiar velocities of the individual stars) is given by 

v r = ARsin2Z, (26) 

where A = (— l/2){Rduj / cLR)q, is Oort's A constant, I is 
the galactic longitude, and for simplicity we have written 
only the expression valid for galactic latitude 6 = (see 
e.g. [25]). 

If the Galaxy is loosing mass, stars become less and 
less bound and acquire radial velocities with respect to 
the Galaxy rest frame. The physics of the effect can be 
understood easily limiting ourselves to a star in a circular 
orbit of radius r in the galactic plane, with a Keplerian 
longitudinal velocity vi given by vf — GM/r; M is an 
effective mass related to the mass at the interior of r. 
In particular the effective mass inside the solar circle is 
obtained from vi ~ 220 km/s and r ~ 8 kpc and is M ~ 
8 x 1O 1O M0. The conservation of the angular momentum 
J = mrvi gives rvi + rvi = 0, or v r /r = —vi/vi, where 
v r = r is the radial velocity in the frame of the Galaxy. 
Using the Keplerian value for vi, this gives 



In the frame of the sun, this gives an additional contri- 
bution to the radial velocity (26), 

v r = ARsin2l + KR, (28) 
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where 



and R is the sun-star distance. In the literature on galac- 
tic dynamics the additional term in cq. (28) is known as 
a K-term. (More generally, one can define a K-term as 
any additional term to be added to the right-hand side of 
eq. (26), independent of the galactic longitude, and with 
an arbitrary dependence on R, e.g. a constant. We see 
that in the case of mass loss one gets a K-term growing 
linearly with R). The difficulty of extracting from the 
data the effect due to the mass loss is that the K term 
receives contributions from many other effects of galactic 
dynamics; for example, the value of K extracted from 
young stars at distances R < 0.6 kpc is strongly influ- 
enced by the kinematic peculiarities of the Gould Belt, 
and therefore the value of K depends on the distance of 
the sample, and it turns out to depend also on the age 
of the stars chosen. To have an idea of the systematics 
involved, we observe that ref. [26], using the Hipparcos 
data, finds that for O.lkpc < R < 0.6 kpc, K (expressed 
in km/(s kpc) ) ranges from 7.1 ± 1.4 if one uses stars 
younger than 30 Myr to —5.4 ± 2.3 for stars with age 
between 60-90 Myr, with an average over all stellar ages 
K = 0.5 ±0.9. Samples including only stars at larger dis- 
tances, and therefore insensitive to the Gould belt, give 
instead negative values of K; in particular, including only 
stars with 0.6 kpc < R < 2 kpc, and averaging over all 
stellar ages, one finds K = —2.9 ± 0.6 [26]. A negative 
value of K might in principle be due to the influence of 
the spiral arm structure, but even taking into account 
spiral arm kinematics, the sample at R > 0.6 kpc gives 
a negative value, K = —(1 — 3) km s -1 kpc -1 [27], and 
the physical origin of the negative sign, representing a 
contraction rather than an expansion, is not really well 
understood. 

To get a bound on the energy loss to GWs, one should 
in principle be able first of all to understand the other 
mechanisms which are at work and which give a neg- 
ative contribution to K, evaluate them precisely, sub- 
tract them, and see if one is left with a positive value 
of K growing linearly with R. This is beyond the ac- 
curacy of our knowledge of the Galaxy. However, in 
any case we certainly do not expect a fine tuned can- 
cellation between some negative contributions to K and 
the positive contribution due to a possible GW emission. 
Let us assume for definiteness that a positive contribu- 
tion to K from GW mass loss is smaller than 20% of 
the absolute value of K. We use the value measured 
at R > 0.6 kpc, to get rid of effects due to the Gould 
Belt, i.e., K = —(1 — 3) km s -1 kpc -1 , and therefore 
we assume that the positive contribution from GWs is 
smaller than O(0.4) kms -1 kpc -1 . From eq. (29), using 
M = 8 x 10 10 M©, this translates into 

-M < O(3O)M /yr. (30) 

Clearly, the precise value of the bound can be modified 



somehow, changing the level of fine tuning that one is 
willing to tolerate, but eq. (30) gives a first order-of- 
magnitude estimate. 



B. Mass loss and outward motion of the LSR 

Rather than looking at the K term, i.e. at the ex- 
pansion/contraction of the stars within a few kpc from 
the sun, one can investigate whether the local standard 
of rest (LSR) has an overall outward radial velocity, 
as suggested by eq. (27). Long ago Kerr [28] indeed 
suggested an outward radial motion of the LSR with 
a velocity wlsr = 7 km/s, as an explanation for the 
lack of axisymmetry of the galactic rotation curve ob- 
tained with 21 cm surveys, and this effect was attributed 
to mass loss due to GWs in ref. [22]. In more recent 
years the experimental determination of ulsr has not 
become much more clear. Blitz and Spergel [29] from 
21 cm line emission find ulsr = +14 km/s (where the 
positive sign means radially outward); results consistent 
with this value have been found from Cephcid kinemat- 
ics [30] while in ref. [31], from a variety of measurements 
(OH/IR stars, globular clusters, high- velocity stars, plan- 
etary nebulae), it is proposed ulsr = — 1 ± 9 km/s. 
From young open clusters [32] one finds a maximum value 
mlsr = 3 ± 2 km/s. Recent work on OH/IR stars [33] 
gives ulsr = 2.7 ± 6.8 km/s. 

However, even if an outward velocity of the LSR of a 
few km/s indeed exists, it cannot be due to mass loss. 
In fact, the observation of the 21 cm absorption line to- 
ward the galactic center [34] shows that the gas along 
the line-of-sight has a mean radial velocity with respect 
to the LSR of -0.23 ± 0.06 km/s. The absorbing mate- 
rial is probably at 1-2 kpc from the galactic center. A 
radial expansion due to mass loss predicts a radial veloc- 
ity v r <~ r, cq. (27), and therefore, if at the sun location 
r ~ 8 kpc this effect has to be responsible for a veloc- 
ity of the LSR, we should see a comparable difference in 
velocity between us and this gas, Av r = (—M/M)Ar. 
The model of Blitz and Spergel based on the existence of 
a rotating triaxial spheroid [29] manages to escape this 
limit because it predicts a dependence of v r on r rather 
flat between 3 and 8 kpc from the galactic center. For a 
mass-loss model instead v r ~ r and we do not have this 
escape route. 

As in the previous section, there can be in general both 
positive and negative contributions from different physi- 
cal mechanisms to the value Av r — —0.23 ± 0.06 km/s, 
and to extract a bound on mass loss to GW (or, for that 
matter, to any mass loss) we require that no fine tuning 
between different contributions takes place. Again we set 
conventionally at 20% the maximum fine tuning that we 
allow, which means that we say that a positive contribu- 
tion from GWs to Av r , if it exists at all, must be smaller 
than 0(0.04) km/s. Setting the distance between us and 
the gas to Ar — 6 kpc, this gives a bound on the mass 
loss to GWs, {-M)/M < 0(0.04) km s -1 /6 kpc; using 
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again M — 8 x 10 10 , this means 

—M < 0(0.5) M /yr . (31) 

Again, we should repeat the caveat that this value de- 
pends somehow on the level of fine tuning that we allow 
but it is clear that to stretch it to, say, 10M Q /yr, we 
must make rather unnatural fine-tuning assumptions. 

C. Upper limits from globular clusters and wide 
binaries 

Upper limits of the same order of magnitudes have 
been found by Poveda and Allen [23] using globular clus- 
ters as probes. The idea is that, if the mass of the Galaxy 
was much bigger in the past, the orbits of globular clus- 
ters would have been much closer to the galactic nucleus, 
and this close interaction with a very massive central nu- 
cleus would have produced the tidal disruption of the 
cluster. Using in particular the globular cluster Omega 
Centauri (NGC5139) as a probe, and assuming a mass 
loss localized in the galactic center, ref. [23] finds that 
the fact that this cluster still exists today implies that 
the Galaxy cannot have lost more than O(l) x 10 10 M Q 
over its lifetime, implying a bound on a steady mass loss 

—M < 0(1) M Q /yr . (32) 

Actually, Omega Centauri is a rather unusual globu- 
lar cluster because its metalicity distribution, mass and 
shape lend support to the idea that it is the remain of 
a dwarf galaxy which has been tidally stripped by our 
Galaxy [35, 36]. In this case the considerations of ref. [23] 
would not apply, since they assume an adiabatic evolu- 
tion of the orbit. For this reason, the analysis of ref. [23] 
has been repeated in ref. [37] for two more globular clus- 
ters, M92 and M5. For M92 one finds basically the same 
limit as for Omega Centauri, while M5 gives a limit larger 
by a factor of two. 

It should be observed, however, that eq. (32) is ob- 
tained assuming a mass loss concentrated in the galactic 
center. Since the motion of the cluster is sensitive only 
to the mass in the interior of the orbit, if the mass loss 
is distributed uniformly over the galactic disk, the limit 
is relaxed by an order of magnitude, and one finds [37] 

-M < O(lO)M /yr, (33) 

from uj-Cen and M92, and —M < 0(20) M Q /yr from M5. 

A limit comes also from the existence of old wide bina- 
ries, since for a very massive galactic nucleus the galactic 
orbits would have been much smaller than at present, 
therefore the density of stars would have been much 
larger and the dissolution time of binaries due to stel- 
lar encounters correspondingly shorter. From a list of 
11 well observed old wide binaries Poveda and Allen [23] 
find a limit on steady mass loss 

-M < O(lO)M /yr. (34) 



D. Comparison with the Explorer /Nautilus data 

Comparing the limits discussed in these sections with 
the mass loss given in eqs. (24) or (25) we can make the 
following considerations. Taking into account the un- 
certainties in the theoretical bounds, we cannot exclude 
that GW emission can generate a steady mass loss at 
the level of eq. (24) or eq. (25). However, such a large 
mass loss can be reconciled with the bounds discussed 
only invoking a certain amount of fine tuning. It should 
also be mentioned that the galactic disk accretes exter- 
nal matter, so this could partially compensate the mass 
loss. However, the accretion is mostly from the edge of 
the disk and in this case it should not influence much our 
considerations. 

We can therefore summarize as follows the points 
which must be addressed in order to advance an inter- 
pretation of the data in terms of a population of sources 
residing in the galactic disk, (i) The estimated merging 
rates of NS-NS, NS-BH and BH-BH binaries are far too 
small, (ii) Exotic objects, for instance related to dark 
matter candidates, are more difficult to rule out on the 
basis of population synthesis models, simply because of 
our ignorance of their dynamics and evolution. However, 
objects of primordial origin, like primordial BHs, should 
reside in the halo and in this case they would not ex- 
plain a correlation between the orientation of the bars 
and the galactic disk, (iii) Independently of its nature, a 
population of sources which produces a steady mass loss 
to GWs, at the level required to explain the data, over 
a timescale comparable to the age of the Galaxy, might 
be difficult to reconcile with the limits from galactic dy- 
namics discussed in this section. Thus, an evolutionary 
scenario able to explain why the mass loss is larger today 
than in the past would probably be needed. 

An answer to the last point might come considering 
sources localized close to the galactic center, rather than 
distributed in the galactic disk. Stellar dynamics in the 
environment of a supermassive black holes can have com- 
pletely different timescales. Extreme examples are given 
for instance by BL Lac objects, whose luminosity can 
change by a factor of 100 over just a few months. The 
recent discoveries of X-ray flaring with a timescale of one 
hour from the galactic center [38], and of rapid IR flar- 
ing [39] point toward the possibility of an active galactic 
nucleus, even if with a very low luminosity compared to 
typical AGN. In any case, for the galactic center, periods 
of enhanced activity on a timescale of, say, millions of 
years, are certainly possible in principle, so it is possible 
that the activity that we are seeing today is larger than 
the average value over the entire age (t ~ 10 10 yr) of the 
Galaxy, and then the limits discussed above on the mass 
loss of the Galaxy would disappear. This hypothesis has 
a clear-cut experimental test, to be checked on future 
data: as we see from fig. 3, a second peak should appear, 
centered at sidereal hour 14. 
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VI. GW BURSTERS 

In this section we explore the possibility that all events 
come from a single source (or at most a few sources), 
which repeatedly emits GW bursts. We will discuss in 
sections VI A- VIC some possible physical mechanisms 
that could produce such a "GW burster". First we dis- 
cuss how this hypothesis might help in the interpretation 
of the data. 

A motivation for this hypothesis is that it allows us to 
overcome the limits on the mass loss discussed in the pre- 
vious section. In fact, if each event comes from a different 
source, it is very difficult to escape the conclusion that 
the average value of the distance r to the source, which 
appears in eq. (7), is at least of the order of the distance 
to the galactic center, r ~ 8 kpc. Indeed, the problem 
is that even within such a distance there is not, as far 
as we know, a sufficiently large population of candidate 
sources. 

If instead we assume that the source is always the same, 
then the detectors will be sensitive to the closest one, 
and therefore r in eq. (7) can be smaller. Of course, if 
the source that we see is "typical" of its population, in 
the sense that its distance from us is of the order of the 
average distance which could be inferred from the num- 
ber of similar galactic objects, we do not gain anything, 
because the factor ~ r 2 in energy which we gain from 
the fact that the source is close to us is compensated 
by the fact that in the Galactic disk there are ~ 1/r 2 
similar sources. The total energy radiated in GWs by 
the Galaxy would then be the same. However, when we 
observe just a single source, statistical considerations do 
not apply, and the source is where it is. In this case 
it therefore makes sense to consider a lucky situation, 
where we happen to be closer than expected to a source, 
or where a transient phenomenon which lasts only for 
a short period happens to be on. If instead each event 
in the EXPLORER/NAUTILUS data is produced by a 
different source, the conclusion that the average distance 
of the sources will be the typical one of the population 
should be inescapable. 

Placing the source close to us opens the possibility 
that it is at work some less powerful mechanism for GW 
production. For a still rather respectable distance of 
500 pc, eq. (7) gives AE lad ~ 4 x lO~ 5 M c 2 for each 
burst while, taking r ~ 100 pc (which is of the order of 
the distances to the closest known neutron stars) gives 
A-EVad ~ 10~ 6 M Q c 2 . From population synthesis models 
it was estimated that the closest NS is at about 5 pc (see 
ref. [40], pag. 12). For a source at such a small distance 
A.E r ad ~4x lO _9 M c 2 . (However, if the source is too 
close, the fact that it is approximately in the direction of 
the galactic plane would have to be explained by chance) . 
Actually, the local density of old NS is reduced by an or- 
der of magnitude by the fact that their birth velocities 
are large [41], so 10 pc is probably a safer estimate. On 
the other hand, the solar vicinity is enriched with young 
NSs (and possibly BHs) which originate in the Gould 



Belt [42]. 

If the source emits O(200) bursts per year it will ra- 
diate 10~ 2 solar masses per year for r ~ 500 pc and 
2 x 10~ 4 M Q /yr for r ~ 100 pc, and therefore it could be 
a transient phenomenon which lasts for a period between 
a few years and a few thousands years. 

As we have seen in sect. Ill, a single source in the 
galactic plane which emits randomly polarized GWs has 
a response function with two peaks, one close to sidereal 
hour 4, and a second at a position which depends on 
the source location. Therefore the signature for such a 
source would be the emergence of a second peak when 
higher statistics will be available. The position of the 
peaks can be shifted if the source emits polarized GWs, 
and in this case the height of the second peak can be 
smaller than the first. 

A process which emits a large burst of GWs should 
be a cataclysmic event which results in the disruption 
of the source, and therefore it should be difficult to find 
mechanisms that produce a GW burster. However, if 
the source is sufficiently close to us we can explain the 
data with relatively small bursts, e.g. 10~ 6 M Q c 2 for r ~ 
100 pc. At this level, it is possible to imagine galactic 
mechanisms that produce a GW burster. In the next 
subsections we will present some possible realizations of 
this idea. 



A. Accreting neutron stars 

Generally speaking, in astrophysics the presence of re- 
peated activity is often related to accretion onto com- 
pact objects. An example is provided by X-ray bursters. 
These are neutron stars which accrete matter from a 
companion. The magnetic field of the NS is not suffi- 
ciently strong to channel the accreting matter toward the 
poles, and therefore the accretion is spherically symmet- 
ric. Each time a layer of about one meter of material is 
accreted (which happens in a time which, depending on 
the particular star, can be between a few hours and a few 
days) a thermonuclear flash takes place, and is observed 
as an X-ray burst. These sources therefore repeatedly 
emit X-ray bursts. 

Stimulated by this example, we ask whether a NS ac- 
creting at a steady rate can undergo periodically some 
structural changes which are accompanied by the emis- 
sion of GWs. We start from the observation that the 
mass-radius relation of a NS depends on the equation of 
state, but it is always such that the larger is the mass, the 
smaller is the radius. Thus, when a NS accretes mate- 
rial, its new equilibrium radius decreases. If the NS were 
a fluid, a continuous accretion of matter would produce 
a continuous decrease of the radius. However, neutron 
stars have a solid crust, about one km thick and with a 
rigidity which, in the inner part, is huge by terrestrial 
standards. Therefore the radius will rather stays con- 
stant until sufficient material has been accreted so that 
the crust can be broken, and the evolution of the radius 



13 



FIG. 14: A schematic drawing of the evolution of R against 
time for an accreting NS. The dashed line is the evolution of 
an idealized fluid NS while the solid line is the evolution of a 
real NS with a solid crust. 



will rather be a sequence of jumps, as shown schemati- 
cally in fig. 14. 

In each jump a certain amount of energy is released. 
Two questions are therefore important for our purposes: 
how much energy is released per jump, and whether this 
energy can be radiated away in GWs. 

Neutron star perturbations can be described in terms 
of quasi normal modes, which can be excited by many 
possible mechanism: accreting material, crust breaking, 
starquakes, onset of phase transitions, coalescence, etc. 
These perturbations can excite in particular the neutron 
star quasi normal modes described by spherical harmon- 
ics with / = 2, which dominates the emission of GWs. 

Furthermore, if the NS is rotating, even the radial os- 
cillation will induce a time varying quadrupole moment, 
and in this case energy will be liberated in GWs, with a 
damping timescale [43] 
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where P ro t is the rotational period. We have taken 14 km 
as a reference value for R, since this is the typical value 
of the radius for a fast rotating NS in a large range of 
masses, 0.3 < M < 1.2M Q , see e.g. sect. 7.1 of rcf. [44]. 

For a NS spinning with P rot of the order of a few ms, 
tqw is smaller than the damping time due to viscosity, 
and the energy liberated will be radiated away mainly in 
GWs (see also refs. [45, 46]). 

The next question is how much energy is liberated in 
each jump. The source of energy is the potential energy U 
of the NS. In order of magnitude, U ~ (3/5)GM 2 /i? and, 
if the crust breaks when a mass AM has been accreted, 
the energy liberated is 
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The two contributions go in the same direction since a 
positive AM induces a negative AR. However, for typ- 
ical equations of states, especially for rotating neutron 
stars, \AR\/R is smaller than AM/M (see e.g. fig. 7.1 
of rcf. [44]), so we just set \AU/U\ ~ 2AM jM. As an 
order of magnitude 
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for a typical NS. Therefore the energy that is released 
when the crust of a NS collapses under the weight of a 
mass AM can be roughly estimated as _E rac i ~ \AU\ ~ 
0.2 AM c 2 . To estimate the value of AM which induces 
the collapse of the crust we consider the pressure P ex- 
erted by the mass AM, 
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and we equate it to the maximum shear stress <7 max that 
can be sustained by the NS crust. The latter has been 
investigated in the context of pulsar glitches [40, 47-51]; 
the maximum shear stress can be written as [49] 



<7 max — M ^max , 
K 



(39) 



where L ~ 1 km is the thickness of the crust, R is the NS 
radius, fi is the shear modulus and # max is the maximum 
strain angle that the crust can sustain without breaking. 
In the lower part of a 1 km thick NS crust, it is estimated 
that fi ~ 10 30 cgs, while in most of the crust region it is of 
order 2 x 10 29 cgs [49]. Concerning 6* max , experience with 
very strong terrestrial materials g 1VGS Pniax ^ 

10- 2 (only 

after hardening and at low temperatures). However, the 
actual value of 8 max is probably lowered by dislocations, 
and could be of order 6< max ~ 10~ 5 to 10~ 3 [51]. There- 
fore, taking L/R ~ 0.1, we get 
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and correspondingly the maximum energy that can be 
liberated in a single starquake is 
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Taking the most favorable values 
#max ~ 10~ 2 we can arrive to AE ra< \ 
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0.5M Q and 



How- 
ever, if we use the more plausible estimate # max ~ 10~ 5 
to 10~ 3 , as well as M ~ 1.4M Q (which in turn gives a 
smaller value of R, R <~ 10—12 km [44]), we obtain much 
smaller values, A£ rad ~ 10~ 12 to 10~ 10 M©c 2 . 

As we discussed above, even placing the source ex- 
tremely close to us, say r <~ 5 pc (corresponding to the es- 
timated distance to the closest NS), we still need a mech- 
anism which radiates at least 4 x 10~ 9 M Q c 2 . So it seems 
that, even with such an extreme choice for r, it is not easy 
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to explain the events seen by EXPLORER/NAUTILUS 
by the breaking of a neutron star crust. Furthermore, 
a very close accreting NS should have been seen as an 
X-ray source. 

Another possibility, however, is that the starquake due 
to the crust breaking, rather than being the source of 
the GWs, might be the trigger for some more important 
structural changes inside the NS. In particular, the core 
of a NS can perform a phase transition from a hadronic to 
a deconfmed quark-gluon phase. In the literature [45, 46] 
it has been considered the possibility that an accreting 
NS acquires a sufficient mass to perform completely the 
transition from a hadronic core to a quark-gluon core. 
The energy gained in the transition, for a NS with a mass 
M ~ 1.5M Q , is of order 0.15M Q c 2 [46], and it has been 
observed that this large energy can excite quasi normal 
modes and be liberated in a GW bursts (with high ef- 
ficiency if the NS rotates sufficiently fast). Particularly 
interesting appear the discontinuity g- modes [52], hav- 
ing typical frequency in the range 0.5-1.4 kHz and con- 
stituting an unique probe for density discontinuities, like 
the ones induced by phase transitions. One can imagine 
that the phase transition does not take place suddenly 
and completely in the whole core. Rather, each time the 
critical mass AM is accreted, a starquake takes place 
and transforms successive layers of the NS core from the 
hadronic to the deconfmed phase. In this case the 0.15 
solar masses will not be released in a single, very large 
bursts, but rather in a series of bursts. One can imagine 
variants of this scenario in which the phase transition in- 
volves strange quark matter, hybrid quark stars [44], or 
a phase of color crystallization [53]. In all cases, unfor- 
tunately, a computation of the energy liberated by the 
phase transition in a NS layer as a consequence of the 
starquake appears a rather difficult task, so this scenario 
remains quite speculative. 



B. Soft Gamma Repeaters 

Soft Gamma Repeaters (SGRs) are X-ray sources with 
a persistent luminosity of order 10 35 — 10 36 erg/s, that 
occasionally emit huge bursts of soft 7-rays, with a 
power up to 10 42 erg/s, for a duration of order 0.1 s. 
There are three known SGRs in the galactic plane 
(SGR1900+14, SGR1806-20, SGR1627-41) and one in 
the direction of the galactic center, SGR1801-23. Fur- 
thermore, SGR0525-66 is located in the Large Magellanic 
Cloud. 

These objects are understood as magnetars [54-57], 
i.e., as magnetically powered neutron stars with huge 
magnetic fields of order 10 14 — 10 15 gauss. The hypothesis 
that SGRs are neutron stars is supported by the hardness 
and luminosity of the bursts, by the periodic modulation 
of the soft tail and by the fact that for at least three of 
them have been found associations with young (t ~ 10 4 
yr) supernova remnants. Experimental evidence for these 
huge values of the magnetic field comes from the observed 



spin-down rate [58]. The magnetar model connects SGRs 
to another class of enigmatic objects, anomalous x-ray 
pulsars (AXPs), which brings the number of known can- 
didates to about twelve [59]. 

Magnetic field lines in magnetars drift through the liq- 
uid interior of the NS, stressing the crust from below 
and generating strong shear strains. For magnetic fields 
stronger than about 10 14 gauss, these stresses are so large 
that they cause the breaking of the 1 km thick NS crust, 
and this elastic energy is suddenly released in a large 
starquake, which generates a burst of soft gamma rays. 

The statistical behavior of SGRs is strikingly similar 
to that of earthquakes: they obey the power low energy 
distribution of Gutenberg and Richter, and they have 
the same waiting time distribution [60]. The bursts ar- 
rive in bunches, when the crust is yielding to the mag- 
netic stresses. For instance SGR1900+14, after decades 
of quiescence, emitted over 50 detected burst during the 
last week of May 1998, and continued to burst into early 
June. It emitted a giant flare on Aug. 27, 1998 (see be- 
low) and overall during nine months in 1998 it emitted 
over 1000 detected bursts. The source SGR1627-41 in- 
stead suddenly showed up emitting about 100 bursts in 
June- July 1998 [59]. 

Occasionally, truly giant flares have been detected. 
One is the March 5, 1979 event from SGR0525-66 in 
the LMC, which fueled interest in SGR astronomy; this 
source emitted overall 16 bursts until May 1983, when 
the activity ceased. No more bursts have been detected 
from this source since. The other giant flare is the event 
on Aug. 27, 1998 from SGR1900+14 [61]. Another very 
bright bursts was emitted by this source on Apr. 18, 2001, 
and several more commons bursts where detected in the 
following weeks, including another large burst on Apr. 28, 
2001 (see [62] and references therein). Giant flares lib- 
erate more than 10 44 erg (i.e. 10~ 10 M Q c 2 ) in gamma 
rays, and have a longer duration, of order 100 s. In the 
magnetar model they are believed to be produced by a 
global large-scale rearrangement of the magnetic field, 
while smaller bursts are produced by local "crustquakes" 
of the NS. 

The emission of GWs from magnetars has been stud- 
ied in refs. [63, 64]. The breaking of the crust produces 
shear waves [65] with a period of the order of the ms, that 
excite non-radial oscillation modes of the NS, damped 
by the production of GWs with a frequency of the or- 
der of the kHz. For "normal" bursts, ref. [63] estimates 
that the total elastic energy that can be released is of 
order 10 45 erg ~ 5.6 x 10~ 10 M Q c 2 , in very good agree- 
ment with our eq. (41) (since, independently of the agent 
which causes the crustquake, either accretion of magnetic 
fields, this is the maximum elastic energy that can be re- 
leased by the breaking of the crust). Ref. [64] performs 
a detailed analysis of the equilibrium configuration of a 
magnetic polytrope and finds that, in giant flares, the to- 
tal energy released by the magnetic field rearrangement 
can be E tot > 10 47 erg ~ 5.6 x 10~ 8 M Q c 2 , and possi- 
bly even as high as E tot > 10 49 erg ~ 5.6 x lO~ 6 M c 2 , 
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for extreme values of the parameters of the model. This 
energy could by radiated in GWs with high efficcncy. 

The distances to observed SGR are highly uncertain, 
but they should all be at r > 5 kpc. In particular the 
distance to SGR1900+14 is estimated at r ~ 5 — 6 kpc. 
Comparing with cq. (7) we see that, if this is correct, 
the energy requirements are not met since, for sources 
at 5 kpc, we need a process which emits more than 
10- 3 M Q c 2 in GWs. 

Nevertheless, magnetars could become very interesting 
candidates for GW production if it were possible to imag- 
ine that one or more of these sources are very close to us. 
Of course, since no very close source has been detected, 
we must provide a mechanism which forbids the observa- 
tion of electromagnetic emission, while leaving the GW 
emission. One generic possibility that comes to mind is 
that the 7-ray emission might be beamed. This cannot 
be excluded a priori in the case of the "normal" flares, 
where the 7-rays originates from localized cracks in the 
NS crust. 

A second possibility is related to the fact that a critical 
value of the magnetic field £?>3x 10 gauss is required to 
suppress the electron cross section on one photon polar- 
ization state below the Thomson value. This decrease in 
the scattering opacity allows the photon luminosities to 
reach the supcr-Eddington values L <~ 10 4 LEdd observed 
in SGR [56, 66]. Neutron stars with a magnetic field 
below this critical value would be much less bright elec- 
tromagnetically, but might still have sufficient tectonic 
activity to produce a significant amount of GWs. 

A distinctive and testable feature of the hypothesis 
that the EXPLORER/NAUTILUS signals come from 
SGRs is a temporal clustering of the events, since each 
source has long periods of quiescence (years or decades) 
until it suddenly enters in a phase of intense activity, with 
bursts arriving in bounces on the timescale of hours to 
1-2 months. From this point of view, it is interesting to 
observe that, out of the eight events at the sidereal hour 
peak in 2001, two came at the distance of one hour from 
each other (see Table 3 of rcf. [1]). If these events cor- 
respond to real GW signals, it is very unlikely that they 
came from two different sources, and they would rather 
suggest a single source with repeated activity. 

Magnetars are believed to stay in their SGR phase in 
the first ~ 10 4 - 10 5 yr of their life. Then, when the 
star cools below a threshold, the dissipation of magnetic 
activity ceases. This is suggested by the theory [57] and 
also fits well with the estimated ages of the supernova 
remnants which are believed to be associated with SGRs. 
When the dissipation of magnetic energy ceases the NS 
enters the so-called "dead magnetar" phase, but the the- 
ory suggests that these stars remain strongly magnetized. 
It is suggested that a large fraction of all NS, say one- 
half, have indeed been active magnetars [59]. In this 
case, the number of NS highly magnetized but presently 
magnetically inactive would be of the same order as "nor- 
mal" NS. Recalling that the distance to the closest NS is 
estimated to be of order 10 pc, the same order of mag- 



nitude estimate would hold for dead magnetars. If these 
objects maintained a starquake activity, possibly related 
to a residual slow diffusion of the interior magnetic field, 
they would certainly be very interesting candidates as 
GW sources. 



C. Strange quark stars 

Another possible realization of a GW burster is pro- 
vided by the r-mode instability in strange stars (see chap- 
ter 12 of ref. [44] for an introduction to strange stars). 
As discussed in ref. [67], in stars made of strange quark 
matter the r-mode instability has a dynamics quite dis- 
tinct from the neutron star case. In particular, in an ac- 
creting strange star the evolution of the GW amplitude 
generated by the r-mode during its first year of evolu- 
tion consists of a repeated series of bursts on a timescales 
from hours to months (see in particular fig. 4 of rcf. [67]). 
The GW amplitude has been estimated to be (eq. (32) of 
ref. [67], with J = 1.635 x 10~ 2 for an n = 1 polytrope) 
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where P rot is the rotation period and a is the r-mode 
amplitude. During the first year of evolution of a young 
strange star the parameter a performs large oscillations 
from very low values, a <~ 10~ 15 , up to values of or- 
der one. These rapid variations therefore result in a 
series of GW bursts. The amplitude in eq. (42), for a 
distance r <~ 1 kpc, is still too small compared to the 
values h <~ 2.5 x 10~ 18 which corresponds to the NAU- 
TILUS /EXPLORER events when the energy deposed in 
the bar is E s = 100 mK, see cq. (8). Still, as stressed in 
ref. [67], the calculation of the GW signal in this early 
phase of the evolution of the strange star is very model 
dependent and there are large theoretical uncertainties. 

The main problem with this source, however, is that 
within the model discussed in ref. [67] this kind of ac- 
tivity can take place only for very young quark stars, 
i.e., in their first year of evolution. Since for the appli- 
cation to the EXPLORER/NAUTILUS data we need a 
galactic, and even relatively close source, this mechanism 
cannot explain the EXPLORER/NAUTILUS result, un- 
less there is a way to rise significantly the duration of the 
activity, as is the case for SGR. 



D. Search for periodicities 

If the events are all due to different sources, there will 
be no time correlation between events. However, if all 
the events are due to a single (or a few) GW burster, 
it makes sense to explore the possibility that there is a 
periodicity in the arrival times of the data. An almost 
periodic behavior could be generated by many different 
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mechanisms. For instance, X-ray bursters show a regu- 
lar behavior since they accrete at a uniform rate (often 
regulated by the fact that it is at the Eddington limit), 
and go off each time that a given critical mass has been 
accreted. So, if detected, such a periodicity would be re- 
ally a "smoking gun" , showing that the events come from 
a single source. We have therefore looked if a fit of the 
form t n = to + nP can reproduce the arrival time of the 
coincidences reported in ref. [1], for sidereal hours 3 and 
4. The analysis is complicated by the fact that, even if 
such a periodicity existed, most of the values t n would 
correspond to times where either the two detectors were 
not simultaneously on, or they were not well oriented 
with respect to the source. Furthermore, in the two bins 
corresponding to sidereal hours 3 and 4, two events are 
expected to be background, and therefore introduce spu- 
rious correlations. Also, one should not necessarily look 
for an exact periodicity but rather for some regular ac- 
tivity, i.e., something of the form t n = t + nP + At n , 
with At n *C P, since in general the mechanism that pro- 
duces the bursts will not be an exact clock, but rather 
will be related to some physical process, like accretion, 
which proceeds more or less at a steady rate. 

Defining f(oj) = J2n ex P{ — i^n}, where the t n are the 
arrival times of the events (and restricting to the 8 events 
at sidereal hours 3 and 4), for a periodic behavior with 
period P the function |/(w)| should have a peak at u> = 
2-ir/P. Of course, a plot of |/(w)| will show very many 
oscillations due to noise, and the point is whether there 
is a peak which emerges clearly above this noise. Within 
the statistics available, we find that this is not the case. 

To understand more quantitatively what this negative 
result means, we consider one of the highest peak of 
|/(w)| (which corresponds to a period P ~ 4.473 days), 
and we verify that with this value of P we are able to re- 
produce the arrival times of 6 out of the 8 events with a 
precision of ±3 hours. We do not assign any positive sig- 
nificance to this result, since |/(w)| displays many other 
peaks of approximately the same intensity. Rather, we 
can use this result to exclude the existence, in the data, of 
very precise periodicities, i.e., writing t n = t +nP+At n , 
we can say that the data exclude a fit of this form, 
with |At„|<l hr. Further statistics will be needed to 
examine the existence of approximate periodicities with 
\At n \>l hr. 

This negative result still gives useful informations, be- 
cause it allows us to exclude mechanisms that would pro- 
duces an extremely regular series of bursts. 



E. Cosmic strings 

A completely different physical scenario is provided by 
loops of cosmic strings. Cosmic strings are topological de- 
fects which appears in various grand unified theories and 
which might be nucleated at a symmetry breaking phase 
transition in the early Universe (for review, see [69]). It 
has long been known that they can be a source of stochas- 



tic background of GWs of cosmological origin [70] . More 
recently, stimulated by a suggestion of ref. [71], it has 
been shown in refs. [72, 73] that cusps in cosmic string 
loops can emit GW bursts whose amplitude, depending 
on the parameters of the model (in particular, the string 
tension and the number of cusps per loops) can be inter- 
esting for GW detection. 

However, in our case an interpretation in terms of cos- 
mic strings suffers of two main problems. The first is 
that, even with the most optimistic values of the pa- 
rameters, the amplitude for the GW burst in the kHz 
region does not exceed h = 1CP 21 , see fig. 1 of ref. [73]. 
While this can be barely detectable with the target sen- 
sitivities of the LIGO and VIRGO interferometers, it is 
much smaller than the amplitude corresponding to the 
EXPLORER/NAUTILUS data which, for a deposited 
energy E„ = 100 mK, is rather /i~2.5x 10~ 18 . 

Second, cosmic string loops are cosmological objects, 
and therefore there is no reason for obtaining a corre- 
lation between the orientation of the bars and the the 
galactic plane. 



F. Summary 

In this section we have examined the possibility that 
all the 0(200) events per year inferred from the EX- 
PLORER/NAUTILUS observations come from a single, 
or at most a few, GW bursters. We have examined dif- 
ferent possible realizations of this idea and, in general, 
we find that the idea is viable only if the source is very 
close (unless the GW energy liberated in each burst is 
much larger than the estimates that we presented). To 
put the source very close, we need a mechanism which 
shuts off the other form of radiation, like X-rays or 7- 
rays, otherwise such close sources would be very bright 
electromagnetically. 

Accretion onto neutron stars does not seem to offer 
this possibility. However, magnetically powered neutron 
stars (magnetars) have a dynamics which depends criti- 
cally on the value of the magnetic field and might offer 
a broader range of possibilities. Further investigation in 
this direction is certainly worthwhile. 

Magnetars also have the rather special property that, 
after long periods of relatively quite life, they suddenly 
become active and emit signals arriving in bounces on a 
timescale of hours to months. Certainly this will be an 
important point to be tested on future data from long 
runs. 



VII. ACOUSTIC DETECTION OF MASSIVE 
PARTICLES 

Beside the issue of the statistical significance of the re- 
sults of the 2001 run, which will be hopefully settled in 
the near future, it is also crucial to understand whether 
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the excitations of the bars could be due to some phe- 
nomenon unrelated to GWs. One important source of 
concern is the fact that particles crossing a bar can lose 
energy by recoil against the nuclei of the material. This 
will produce a warming up and a local thermal expan- 
sion, which in turns originates mechanical vibrations in 
the bar. The vibrational energy E s deposed in the fun- 
damental mode of a cylindrical bar by this mechanism 
has been studied by several authors [74-80] and is given 
by the formula 
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where L is the bar length, R is the bar radius, lo the 
length of the particle's track inside the bar, zq the dis- 
tance of the track midpoint from one end of the bar, 9 
the angle between the particle track and the axis of the 
bar, dE / dx the energy loss per unit length of the particle 
in the bar, p the bar density, v the sound velocity of the 
material and jc the Griineisen coefficient of the material, 
which depends on the ratio of the thermal expansion co- 
efficient to the specific heat (p/Q ~ 1.6 for Aluminum at 
low temperatures). 

This formula has been well verified experimentally 
for cosmic rays, when the resonant bar is in a non- 
superconducting state (which is the mode in which the 
bar is operated in the 2001 run that we are discussing), 
while it appears to underestimate E s when the bar is in 
the superconductive state [81]. 

Cosmic rays are detected by two detectors placed above 
and below each resonant bar, and the corresponding 
events are eliminated from the list of GW candidate 
events. However a massive, electrically neutral, long- 
lived particle could avoid the cosmic ray veto and depose 
energy in the fundamental vibrational mode of the bar, 
according to eq. (43). In this section we therefore explore 
whether it is possible for massive particles to produce a 
signal compatible with the EXPLORER/NAUTILUS ob- 
servations. 

If the existence of a correlation between the energies 
deposed in the two bars will be confirmed by further data, 
it will be extremely difficult to see how such a result could 
be originated by particles whizzing around, such that oc- 
casionally one particle hits one bar while a second particle 
hits the other bar. Observe in particular, from table 3 of 
rcf. [1], that for each of the eight events at sidereal hours 
3 and 4, the energy £?Ex P i deposed in EXPLORER and 
the energy E^ au deposed in NAUTILUS always differ by 
less than 20% and in some cases by about 5%. From 
one event to the other, instead, the energy can change 
by as much as a factor of 4. Similarly, it is difficult to see 
how a sidereal time modulation could emerge from the 
interaction with a random flux of particles. 

A dependence on sidereal time could instead be under- 
stood if the two bars are detecting particles which are 



ejected in astrophysical explosions. Even in this case, it 
is not easy to explain the existence of a strong correla- 
tion in energy. One should make the assumption that the 
particles are produced with an energy spectrum which is 
strongly peaked, so that the typical spread in energy is 
small, with AE/E <~ 5 - 10%. The position of the peak 
in the energy spectrum might depend on the conditions 
in the star, so it could change from event to event (as it 
is observed in the data), still maintaining a correlation 
between the energies detected in the two bars. Such a 
scenario does not seem very plausible; nevertheless, let 
us assume it for the moment and see if we can find more 
quantitative arguments in order to rule it out. 

First of all observe that in this scenario the particles 
come from astronomical distances and are massive, so any 
spread Aw in their velocities will produce a corresponding 
spread At in the arrival times. Since the bursts seen by 
the bars are concentrated in an interval At < 0.1 s, the 
particles must be highly relativistic, so that their spread 
Av is sufficiently small. In terms of 7 = (1 — w 2 /c 2 ) -1 / 2 
we have, for v close to one, 
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Since A7/7 = AE/E and we are assuming that 
AE/E<1 in order to account for the energy correlation, 
we have Av/v ~ I/7 2 and therefore At/t ~ I/7 2 , where 
At is the duration of the bursts, and t ~ r/c is the time 
taken by the particle to arrive from a distance r. Requir- 
ing At < 0.1 s gives 
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Therefore, for astronomical distances, we must have 
highly relativistic particles. 

Observe that eq. (43) has a strong angular dependence 
which is due to the fact that a particle which crosses the 
bar along its longitudinal axis has a longer path inside 
the material, and therefore a larger energy loss, com- 
pared to a particle which goes through the bar transver- 
sally. When a particle enters at an angle 9 = (i.e. in 
the longitudinal direction) we have lo = L,zq = L/2, 
and the term [. . .] 2 in eq. (43) becomes L 2 /(irR) 2 . On 
the other hand, for a particle entering perpendicularly 
to the bar we have cos 6* — > and Iq = 2R so, in 
eq. (43), [. . .] 2 — ► sin 2 (7rz /L) — ► 1/2 (averaging the sin 
squared over zq). Therefore the energy deposed by parti- 
cles impinging on the bar longitudinally is larger than for 
particles impinging transversally by a factor 2L 2 / (irR) 2 , 
which is of order 20 with the values L ~ 3 m, R ~ 0.3 m 
of the bars. The dependence on 9 of the response func- 
tion, shown in fig 15, is quite peaked, since the maximum 
angle for which the particles traverse the bar in the lon- 
gitudinal direction is 9 — 2R/L ~ 0.2. Therefore, even 
as an acoustic detector of particles, the bar is sensitive 
to the orientation with respect to the source. 



18 



0.8 



0.6 



0.4 



0.2 



0.2 0.4 



0.6 



0.8 

e 



1.2 



1.4 



FIG. 15: The sensitivity curve of the bar for the acustic detec- 
tion of massive particle, as a function of the incoming angle 
9 (solid line). To appreciate the fact that the response func- 
tion is quite narrowly peaked, we also plot for comparison the 
function cos 2 6 (dashed line). The cusp at 6 = 2R/L ~ 0.2 
reflects the fact that this is the maximum angle for which the 
particles traverse the bar longitudinally. We averaged over all 
possible impact points on the surface of the bar. 



It is useful to separate the discussion into two parts, 
corresponding to particles traversing the bar along the 
longitudinal or the transverse direction. 



A. Longitudinal trajectories 

Since the excess of coincidences is seen when the bar 
is orthogonal to the galactic plane, particles traversing 
the bar longitudinally would come from the directions of 
the galactic poles. So we should ask what lies in these 
directions and an answer is that, very close to the North 
Galactic Pole, there is the Virgo cluster of galaxies. This 
result is quite interesting because, with the 2500 galaxies 
present in the Virgo cluster (which include giant galaxies) 
a rate of order 200 events per year could be explained 
by explosions taking place at a rate of a few events per 
century per galaxy, quite compatible, for instance, with 
the rate of supernova explosions. 

In fig. 16 we show the sensitivity curve of the bar as 
an acoustic particle detector, i.e., we show the energy de- 
posed by a particle in the bar as a function of sidereal 
time, for a fixed energy of the incoming particle and a 
given source location. The sensitivity to sources located 
in the Virgo cluster (dashed line) is peaked at sidereal 
hour 3.9 (using as always the average sideral time be- 
tween EXPLORER and NAUTILUS). For comparison, 
we also show the sensitivity to sources in the galactic 
center (solid line), which has a maximum at sidereal time 
21. In both cases, the peaks are quite narrow. The sen- 
sitivity curves have been normalized so that the value at 
the peak is one. 

Inserting the numerical values in eq. (43) we find that 
a signal E s ~ 100 mK can be obtained from a particle 
which crosses the bar longitudinally, with an energy loss 




t 



FIG. 16: The response function of the bar for the acustic 
detection of massive particle, for a source located at the center 
of the Virgo cluster (dashed line), and for a source in the 
galactic center (solid line), as a function of sidereal time t. 
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dE 
dx 

or, in a generic material, 
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Traversing the bar longitudinally (L = 3 m) the particle 
will then lose approximately 6 x 10 3 GeV, and a simi- 
lar amount will be lost traversing the atmosphere (since 
Pai x(3m)~ pAir x (10km)); therefore the particle must 
have had an initial energy of at least E ~ 10 4 GcV. This 
energy arrives on a surface ttR 2 ~ 0.3 m 2 . The accep- 
tance of the detector for particles traversing in the lon- 
gitudinal direction is Al = irR 2 8Q, where <50 is the solid 
angle subtended by the particles. For particles travers- 
ing the bar in the longitudinal direction the solid angle is 
quite small, Stt ~ (2R/L) 2 sr and A L ~ 1.1 x 10~ 2 m 2 sr. 
Then the total energy radiated in the burst must have 
been 



E t , 



E 



Anr z 



10 



-li 



E 



.) 



10 4 GeV7 Vlkpc 



M Q c 2 . 
(48) 

In order to obtain an astrophysically realistic value we 
require that E to t < 10~ 2 M©c 2 . This gives 



1< 



E 



10 4 GeV ' 



< 



32 Mpc y 



(49) 



where the lower bound comes from the fact that and ini- 
tial particle energy E> 10 4 GeV is needed in order to de- 
pose 100 mK in the bar. Eq. (49) also implies a maximum 
value for the distance to the source, r<32 Mpc; this dis- 
tance include the Virgo cluster, which is approximately 
at r <~ 16 Mpc. Writing E — ■ym and using eq. (45) we 
also find from eq. (49) a limit on the particle mass, 

5/2 

m < 10 7 f ) GeV . (50) 
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yr x 1.1 x 10~ 2 m 2 sr 



For a source in the Virgo cluster at r = 16 Mpc this gives 
to < 0.3 MeV. 

Our final information on these hypothetical particles is 
on their flux. The event rate, as we have seen, is 0(200) 
events per year. Using the acceptance At found above 
the flux should be 

200 „ „„ o o i i 

~ 6 x l(r 8 cm _2 s _1 sr _1 . 

( 51 ) 

Summarizing the above discussion, one might formulate 
the hypothesis that the EXPLORER/NAUTILUS data 
could be due to electrically neutral massive particles com- 
ing from explosions in the Virgo cluster. In this case the 
sidereal time dependence could be explained, and one 
could as well find a sufficient number of sources to jus- 
tify the observed rate. The candidate particle should 
have 7 > 10 8 (from eq. (45) with r = 16 Mpc), a mass 
to < 0.3 MeV, an energy 10 4 GeV < E < 4 x 10 4 GeV, 
a stopping power (1/ ' p)dE/dx = 7GeVcm 2 /gr and a flux 
F ~ 6 x 10~ 8 cm^s^sr- 1 . 

The existence of a particle with these properties is 
however excluded. WIMPS are obviously out of ques- 
tion: with their typical energy loss on nucleons of or- 
der 100 keV and weak cross sections, WIMPS have 
(\/p)dE/dx <~ 10~ 27 GeVcm 2 /gr. More in general, 
if a particle with a stopping power (1/ p)dE/dx = 
7GeVcm 2 /gr and a flux F ~ 6 x 10~ 8 cm _2 s _1 sr _1 ex- 
isted, it would have been detected by experiments on 
cosmic rays at altitude, which for these stopping pow- 
ers are sensitive to fluxes many orders of magnitudes 
below our required value. The precise value of the ex- 
perimental upper limit on the flux depends on the spe- 
cific properties of the particle which is being searched. 
However, to get an idea of the order of magnitudes, we 
may observe that cosmic ray experiments have put lim- 
its on the fluxes of hypothetical objects like nuclearites 
(bound states of u, d, s quarks surrounded by an electron 
cloud) and Q-balls. These exotic objects are extremely 
massive and are therefore excluded as our candidate par- 
ticle. However their energy losses in matter, which is 
the main parameter for detection, are comparable to the 
value that we need. The upper limit on the flux of these 
particles are of order 10" 



experiments and 10 



46 cm 2 s 



sr from underground 



detectors like MACRO (see ref. [82] for review) 

Furthermore, a strongly interacting particle with a 
mass below the MeV would spoil primordial nucleosyn- 
thesis. 



B. Transverse trajectories 

Recall from the discussion of sect. Ill A that the sensi- 
tivity curve determines the distribution of the number of 
events versus sidereal time only when the typical SNR 
is low. In the opposite limit where almost all events 
are above threshold, an event will be detected indepen- 
dently of its arrival direction. In this case the bar will 



detect more events when the particle flux is transverse 
to the bar axis simply because the geometric acceptance 
for transverse particles is much larger than for longitu- 
dinal particles. In fact, the geometric cross section is 
now S = 2RL ~ 1.8 m 2 and the solid angle is of order 
47r (more precise estimates are not necessary for our or- 
der of magnitude calculation). Therefore the acceptance 
for particles arriving transversally is At — An(2RL) ~ 
22.6 m 2 sr, to be compared with the acceptance for par- 
ticles arriving longitudinally, At — 1.1 x 10~ 2 m 2 sr, as 
found in the previous subsection. 

In the case of sources with high SNR in the galactic 
plane the dependence of the number of events on 9 would 
then be proportional to sin#, because Lsm.9 is the pro- 
jection of the bar length in the direction orthogonal to the 
galactic plane. In a two-detector correlation the number 
of coincidences will be proportional to sin 2 9; this func- 
tion has a maximum around sidereal hour 4, so the loca- 
tion of the peak can be reproduced, but the peak in this 
case is rather broad, so this scenario is is not favored by 
the present data, which rather suggest the presence of a 
relatively narrow peak, concentrated in two hours. 

The analysis of the possible particle physics candidates 
proceeds similarly to the previous sections, with two main 
differences: first, because of the angular dependence of 
eq. (43), if the particle goes through the bar transversally, 
in order to depose 100 mK in the bar we need a value 
of dE/dx larger by a factor \/20 ~ 4.5 compared to the 
longitudinal case, i.e. 



dE nn GeV 
-90 

dx cm 



(52) 



in Aluminum. For a transverse length 2R = 60 cm this 
gives a loss of 5.4 x 10 3 GeV in the bar, very similar 
to the longitudinal case, plus about the same loss in the 
atmosphere, so again the initial particle energy should 
have been at least 10 4 GeV. 

A more important modifications is due to the fact that 
the acceptance At of the detector is now much bigger, 
and eq. (48) becomes 



Etnt — E 
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5x10" 
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M P) c 2 



10 4 GcV; Vlkpc; 

(53) 

Correspondingly, eq. (49) is replaced by 



1 < 



E 



10 4 GeV 



< 2 x 10 



12 



1 kpc 



(54) 



where, beside using the new acceptance, we also used 
a typical galactic distance as a reference value for r. 
Eq. (50) becomes 



5/2 



. i GeV. 



(55) 



The flux would instead be 
200 



yr x 22.6 m 2 sr 



~ 3 x 10" 11 cm"^"^:" 1 . (56) 
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A flux of this order of magnitudes still remains very dif- 
ficult to reconcile with cosmic ray experiments (even if 
some uncertainty can appear, because the experimental 
limits depend on the detailed interaction properties of 
the particle in question). 

Furthermore, despite the fact that the constraint on 
the mass given in eq. (55) is now much less tight, still 
there is no plausible particle physics candidate. In par- 
ticular, hypothetical neutrons coming from an astrophys- 
ical source would be stopped in the atmosphere in about 
1 km; all neutrons detected at see level are produced in 
the collision of primaries with atmospheric nuclei, so they 
cannot account for a correlated signal in the two bars. 
Exotic objects like nuclearites and Q-balls are still ex- 
cluded by the experimental limits on the flux, which are 
3 orders of magnitude tighter than the value required by 
eq. (56). More in general, any object coming from outer 
space, with an energy loss of the order of the value given 
in eq. (52) and a flux at see level as in eq. (56) should 
have been seen in cosmic ray experiments at altitude. 

In conclusion, it seems that we can safely exclude that 
the coincident events are caused by the passage of mas- 
sive particles. 



VIII. CONCLUSION 

In the EXPLORER/NAUTILUS 2001 run there is a 
potentially interesting effect, i.e., an excess of coinci- 
dences when the two detectors are perpendicular to the 
galactic plane, with a correlation in energy between the 
events in the two bars. Therefore, we think that it is ap- 
propriate to start investigating the possibility that galac- 
tic burst sources could explain these observations. It is 
evident from our discussion that there is no obvious phys- 
ical explanation for these data. We have examined a 
number of different possibilities, and discussed the dif- 
ficulties of each. If confirmed, a crucial clue for un- 
derstanding this effect will be provided by the sidereal 
time distribution of events, which, with sufficiently large 
statistics, will be able to discriminate between a popula- 
tion of sources distributed in the galactic disk, a popula- 
tion of sources located at the galactic centers or a single 
close source. 

Another important clue might be a time clustering 
of the events, with recurring bursts on the timescale of 
hours, and periods of highly enhanced activity on the 
scale of the month. This would suggest a kind of activity 
similar to what is observed in Soft Gamma Repeaters. 

New data from GW detectors (the large interferome- 
ters LIGO, VIRGO, GEO and TAMA and the bars AL- 
LEGRO, AURIGA, EXPLORER and NAUTILUS) will 
soon be available, and hopefully the situation will be clar- 
ified in the near future. 




galactic plane 



FIG. 17: The geometric setting discussed in appendix A. 

APPENDIX A: COMPUTATION OF 0(i),ip(t) 

We have seen in eq. (9) that the response function of 
the detector depends on the angle 9 and ip. In turn, these 
angles depend on sidereal time because of the rotation of 
the Earth. The dependence of 9 on t has been given 
in eq. (10). In this appendix we perform the explicit 
computation, both for 9{t) and for ip(t). 

For definitcness, we consider a source in the galactic 
plane, but the computation is easily adapted to the most 
general case. The geometric setting is shown in fig. 17. In 
the galactic plane, we consider a source S and the Earth 
(marked by E). On the source, we build a reference frame 
with the x' axis pointing toward the Earth and the z' 
axis in the direction of the Galactic North Pole (GNP). 
Therefore the y' axis lies in the Galactic plane. The GW 
propagates along the x' axis and therefore its transverse 
plane is (y 1 , z'), and the + and x polarizations are defined 
with respect to the y' and z' axes. 

On the Earth we build a reference frame with the x axis 
pointing toward the source, the z axis toward the GNP 
and therefore the y axis in the galactic plane. We denote 
by d the unit vector in the direction of the longitudinal 
axis of the bar. The angle between d and the x axis is 9. 
The angle ip instead describes how the bar is rotated in 
the (y, z) plane, i.e., with respect to the axes which are 
parallel to the (y 1 , z') axes used to define the polarizations 
+ and x . Therefore the unit vector d can be written, in 
the (x, y, z) frame, as 

d = x cos 9 + y sin 9 cos ip + z sin 9 sin ip . (Al) 

On the other hand, the direction d of the bar is given 
by its location on Earth (expressed by the latitude I and 
longitude L), its azimuthal angle a and the local sidereal 
time t. At the detector location, we define a new coor- 
dinate system with axes (E, Nl, Z), where E is the east 
direction, Nl the local north (both vectors are therefore 
in the plan tangent to the Earth surface) and Z points 
toward the zenith, so —Z points toward the Earth center. 
In this frame the bar direction is given by 

d = E sin a + N L cos a . (A2) 

The position of a source in the sky, with respect to the 
Earth, is given in equatorial coordinates, i.e., in the basis 
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FIG. 18: The rotation between the coordinates systems 
(E,N L ,Z) and (E,N,E±). The unit vector E is perpen- 
dicular to the plane of the figure, pointing downward. 



GNP 




FIG. 19: The vector v = N x z, where z point in the direction 
of the Galactic North Pole. The vector v is in the (7,7±) 
plane. 



(7, 7_L,-/V), where N points to the celestial north pole, 7 
to the vernal point and A /± — N x 7. In this basis a unit 
vector n pointing toward a source is expressed in terms 
of the angles (a, 5) as 

n = 7 cos a cos <5 + j± sin a cos <5 + TV sin <5 . (A3) 

In order to express d in equatorial coordinates, we first 
pass from the basis (E,N^,Z) to the basis (E, N, E±) 
where E± = E x N . This is done with a rotation by 
an angle I around the direction E, see fig x 18. Next we 
perform a rotation by an angle fit around N. This brings 
us to the basis ( , y±,N, 7). (If instead of working with 
local sidereal time t we want to use Greenwich sidereal 
time, we must rather rotate by an angle f^Green. + L). 
Therefore in the basis ( A f±,N, 7) the components of d are 
given by 



(cos fit sin fit ' 
1 

— sin fit cos fit , 

or, explicitly 





cos I sin I 
- sin I cos I 




d = 7 (— sin a sin fit — cos a sin? cos fit) + (A5) 
+ 7_l (sin a cos fit — cos a sin I sin fit) + N cos a sin I . 

Using eqs. (A3) and (A5) we can now compute cos#(i) 
from cos9(t) — n ■ d, and we get eq. (10). 

In order to compute ip(t) for a source in the galactic 
plane we rather need to transform eq. (A5), which gives 
the ^-dependence of d, into the (x, y, z) frame, and then 
we can read ip by comparing with eq. (Al). Therefore we 
must find the rotation that brings the frame (7, 7_i_,-/V) 
into the frame (x,y,z). This can be done in steps, in- 
troducing first a vector v = N x z. The unit vector 
v = v/ || v || is given by 



v = —7 sin cvqnp + 7± cos «gnp 



(A6) 



where «gnp is the right ascension of the Galactic North 
Pole (recall that z points in the direction of the GNP). 



The vector v lies in the plane (7,7.1), sec fig. 19, and 
makes an angle f}\ with the 7 axis, with 



(5\ = arccos(7 • v) — arccos(— sinaQNp) • 



(A7) 



Then we can pass from the frame (7,7j.,iV) to the 
frame (v,v±,N), where v± = N x v, performing a rota- 
tion around N by an angle f3\. Next we rotate the N 
axis into the z axis, performing a rotation around v by 
an angle 



02 — arccos(A • z) — arccos(sin Sgnp) ■ 



(A8) 



This rotation leads into the basis (v, v'±,z), where v and 
v'± = z x v are both in the galactic plane. 

Finally, we perform a rotation in the galactic plane in 
order to bring the first axis in the direction of the desired 
source. This is accomplished by a rotation around z by 
an angle (3%. If n is the direction of the source and (a, 6) 
its right ascension and declination, then 

cos Ps — n ■ v = — sin cvgnp cos a cos S + 

+ cos «gnp sin a cos S . (A9) 

In conclusion, the components of d in the frame (x, y, z) 
are given by 



d x \ 


/C3 


dy 


= [ss 


dj 


\ 








(A10) 



where Cj = cos/?i, Sj = sin/Jj, (i = 1, 2, 3) and rf 7 , d 1± ,dN 
are the component in the (7, ^f±,N) frame. Comparing 
with eq. (Al), we then find tp(t) from tan^ = d z /d y . 
Figs. 9-11 are obtained from these expressions, shifting 
sidereal time in order to obtain the mean sidereal time 
between the NAUTILUS and EXPLORER locations. 
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FIG. 20: The geometry for microlensing discussed in the text. 
S,L,0 denote the source, lens and observer, respectively. 



APPENDIX B: MICROLENSING OF GWS 

A magnification effect could help in explaining the re- 
sults, and we have investigated whether the microlensing 
of GWs can play a role in our problem. We find that the 
answer is clearly negative, but still we consider interest- 
ing to resume in this appendix the reasons. 

Gravitational waves, just like electromagnetic waves, 
can be lensed by a large mass situated between the source 
and the observer. The standard calculation of microlens- 
ing is performed using geometrical optics, and in this 
approximation the amplification factor A in the energy 
density is (see e.g. ref. [83]) 



A = 



1 



iVu 2 + 4 u 



(Bl) 



where u = (3/9e, (3 is the angle of the source with re- 
spect to the observer- lens axis (see fig. 20) and 9e is the 
Einstein angle, 



2R S 



Dol(D S l + Dol) 



1/2 



(B2) 



Here Rs is the Schwarzschild radius of the lens and D$l 
and Dol are defined in fig. 20. The second equality in 
eq. (Bl) holds when u -C 1, i.e. when the source, lens 
and observer are well aligned. 

There is however a crucial difference between the am- 
plification of electromagnetic waves and of GWs. The 
geometric optics approximation holds when the reduced 
wavelength of the wave, X/(2tt), is much smaller than 
the typical curvature radius of the spacetimc, which in 
this case is given by the Schwarzschild radius of the 
lens, Rs- When \/(2ir) becomes of order Rs diffraction 
effects become important and the magnification disap- 
pears [84-87]. When the lens is a stellar mass object, its 
Schwarzschild radius is of the order of a few kms. For 
visible light, therefore, the condition X/(2n) <C Rs is 
very well satisfied. On the contrary, the GWs searched 
at resonant bars have a frequency / — 1 kHz, and there- 
fore A/(2-7t) <~ 50km 3> Rs- This means that stellar 
mass objects do not amplify GWs of a detectable fre- 
quency To obtain some amplification, we need a lens 
with a Schwarzschild radius of at least hundreds of kms 
so, within the Galaxy, the only possibility is given by 
supermassive BHs. In this case A/(27r) <C Rs, and the 



maximum amplification is given by [85] 



Ad 



4ir z 



A 



(B3) 



In other words eq. (Bl) saturates, because of diffraction 
effects, at a value of the angle (3 = (3* such that 



(B4) 



The supermassive BH at the center of the Galaxy has 
M~3x 10 6 Mq, and therefore it could provide a maxi- 
mum amplification factor 



^max — 1 X 10 . 



(B5) 



Unfortunately, this factor is of no help for our purposes. 
In fact, using eq. (B2) for 9e, setting Dol — 8 kpc, 
which is the approximate distance to the galactic center, 
and using for simplicity the approximation D$l *C Dol 
(the result below can change at most by a factor of two 
without this approximation, when the source is still in 
the Galaxy), we find that such an amplification is reached 
only if the source, lens and observer are aligned within 
an angle (3 < (3* with 
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and correspondingly the distance ry between the source 
and the observer-lens axis (see fig. 20) must be smaller 
than 



77* ~ 3 x 10" 



D 



SL 



pc 



1/2 



R(7) . 



(B7) 



This is a ridiculously small distance, and the chances of 
finding a source (not to mention one hundred sources 
every year!), which emits a GW burst just when it is so 
precisely aligned between us and the central BH are zero. 
Smaller magnification factors could be obtained within 
a larger region, but, besides the fact such events still 
remain extremely unlikely, smaller magnification factors 
are also of limited help because in the case of microlensing 
the source must be behind the central BH, and therefore 
at a distance from us r > 8 kpc. Clearly, to alleviate 
the energy requirement in eq. (7) by, say, a factor of 
100 it is much easier to assume that the source is at a 
smaller value of r, say r<l kpc, rather than hiding it in 
an extremely narrow cone beyond the central BH. 

Another possibility that one can consider is that there 
is a source of continuous GWs, rather than of bursts, that 
orbits the central BH, in such a way that at a certain 
moment the source is well aligned for microlensing, and 
its emission is suddenly amplified. The idea however does 
not work since, even assuming an amplification factor 
10 6 , in order to produce a burst of an apparent energy 
10~ 2 Mqc 2 , released in a period of less than 0.1 s (which 
is an upper bound on the duration of the bursts seen by 
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the bars) the source should have a continuous emission 
at a rate 10 _8 M Q c 2 in a time < 0.1 s, corresponding 
to a steady rate > 3M Q c 2 /yr. No continuous known 
source can emit at such a huge rate, not to mention the 
implausibility of the hypothesis that the location of the 
orbit is so well fine-tuned for performing microlensing. 

In conclusion, microlensing of GWs does not appear to 
play a role in our problem. 
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